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We formulate the Dynamical Mean Field Approximation equations for the double-exchange 
system with quenched disorder for arbitrary relation between Hund exchange coupling and electron 
band width. Close to the ferromagnetic-paramagnetic transition point the DMFA equations can be 
reduced to the ordinary mean field equation of Curie-Weiss type. We solve the equation to find the 
transition temperature and present the magnetic phase diagram of the system. 
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I. INTRODUCTION 

The Dynamical Mean Field Approximation (DMFA) 
(see and references therein) is a widely used formal- 
ism for many body systems. In many cases it well takes 
into account non-trivial effects of strong correlations and 
thermal disorder. In particular, for the double-exchange 
(DE) model [^|J|], where the exchange interaction be- 
tween the core spins is mediated by mobile carriers, appli- 
cation of the DMFA to the description of ferromagnetic 
transition (see Ref. |Q and references therein) resolved a 
long standing problem of the energy scale for the Curie 
temperature T c of manganites. However, this application 
was restricted to the case where only thermal magnetic 
disorder is present in the model. Thus, chemical disorder 
introduced by doping impurities, which is generic for the 
manganites and many other DE systems, was ignored. 

It was shown however, that the concurrent action of 
the static and magnetic disorder is crucial for the the de- 
scription of the quasi-particle and transport properties of 
the DE systems Therefore the inclusion of static 

disorder into the DMFA scheme is actual. The present 
note is devoted to this task. 



II. HAMILTONIAN AND DMFA EQUATIONS 

Consider the DE model with random on-site energies. 
The Hamiltonian of the model is 

H ^ ^ ^n — n' clia^na "F ^ ^ V n C^ not C n ' a 
nn' a na 

- m„ • (J a ^ na Cn&, (1) 
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where t n _ n i is the electron hopping, V n is the random on- 
site energy, J is the effective exchange coupling between 
a core spin and a conduction electron, <x is the vector of 
the Pauli matrices, and a, (5 are spin indices. We express 



the localized (classical) spin by m„ = (m n x ,m n y ,m n z ) 
with the normalization |m| 2 = 1. 

In a single electron representation the Hamiltonian can 
be presented as 



H nn > = H n _ n , + (V n — Jra„ • a) 5 nn >; 



(2) 



the first is translationaly invariant, the second describes 
quenched disorder, and the third - annealed disorder. 

The DMFA, as applied to the problem under consider- 
ation, is based on two assumptions. The first assumption 
is that the averaged, with respect to random orientation 
of localized spins and random on-site energy V , locator 



where 



G\ oc (z) = (G nn (z) 



G{z) = {z-H)-\ 



(3) 



(4) 



can be expressed through the the local self-energy £ by 
the equation 



G 



ioc(z) = go (z- t(z)j , 



where 



«.w = ^E( z - fl £)" 1 



(5) 



(6) 



is the bare (in the absence of the disorder and exchange 
interaction) locator. Thus introduced self-energy satisfies 
equation 



Gloc(z) = v _ . 



1 



(7) 



The system of equations (^|) and (j^) is very much similar 
to the well known CPA equations (see B and references 
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therein), as generalized to the case when the quantities 
G, £ and g are 2x2 matrices in spin space . The 
system of equations however, is not yet closed. The av- 
eraging with respect to annealed disorder is principally 
different from the averaging with respect to quenched 
disorder. 

The second assumption of the DMFA is the prescrip- 
tion for the determining, in our case, the probability of 
a spin configuration self-consistently with the solutions 
of the Eqs. (j|) and (0). To formulate the DMFA equa- 
tion for this probability, taking into account both kinds 
of the disorder, let us start from the general formula for 
the partition function 



^xp ^-Tr^logG(z s )^ ndm„, 



(8) 



where z s = ioj s + /i; lo s is the Matsubara frequency and 
fi is the chemical potential. The averaging over {m„} is 
given by 



linear with respect to magnetization M approximation 
Eq. (|l^) takes the form 

P Vn (m) oc exp (-/3I Vn M ■ m) . (14) 

Non-trivial solution of the MF equation 

M = y (P V7i (m)} v mdm. (15) 
can exist only for T < T c , where T c — ^ {Iv n )v 

III. T c FOR THE SEMI-CIRCULAR DOS 

The problem of finding T c is reduced to calculation 
of gv and hy- For simplicity consider the semi-circular 
(SC) bare density of states (DOS) N (e), the bandwidth 
being 2W. Then 



9o(z) 



N (e)de _ 2 
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For this case 



£ = z - 2wG\ oc - G lo ]., 



(16) 



(17) 



All observables, in particular thermodynamic potential 
f2, should additionally be averaged over the realizations 
of the quenched disorder; in particular 



n 



(10) 



where w = W 2 /8. Thus from Eqs. ([}]) and (m) we obtain 
a single equation for G\ oc 



(18) 



The DMFA approximates the multi-spin probability 
TrlogG) as a product of one-site probabil- 
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1 exp (- 



ities in such a way, that 
5Q 



SG \oc 



0. 



(11) 



The result for the one-site probability reads (for details 
of the calculation see Ref. Jl0|): 



Gi oc (z) = ( s ) 

\z- 2wG\ oc {z) - V n + Jm • a j m y 

andEq. (|l|) can be presented as 

z s - 2wG\ oc {z s ) 

-V n + Jm-a} e iWs0 +. 
In linear with respect to M approximation 
Gioc = gl - hJM ■ a, 



AO ul i , - ij^logdet 



(19) 



(20) 



iV„(m) oc exp [-/?ASl m) yJ 



(12) 



where g is locator in paramagnetic phase, given by the 
equation 



where 



m ,y„ = — ^THog [l + G loc (z s ) 



1 
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(jm • a - V n + t(z s ) 



Jw 3 + 



(13) 



z - 2wg -Vn-J/ V 
1 



is the change of the thermodynamic potential of the elec- 
tron gas described by the Green's function Gi oc due to 
interaction with a single impurity [O 12 1. 

The right hand side of Eq. ( |l2| ) , is a complicated non- 
linear functional of Py (m). However, if we are interested 
only in the transition temperature T c , the problem can 
be reduced to a traditional mean field (MF) equation. In 



, z - 2wg - V n + J i VJ 
and the quantity h is given by the formula 



1 



where 



Ay„(z s ) = 



1 



\z s ~ 2wg(z s ) - V n ] - J 2 



(21) 



(22) 



(23) 
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Expanding Eq. ( jl9| ) we obtain the effective exchange in- 
tegral is 



_ AJ 2 w \ , , . . . . 



(24) 



If we transform the sum over the imaginary Matsubara 
frequencies in the right-hand side of Eq. (p4|) to integral 
over real energies E, we obtain for the T c 

4 T 2 w r°° 

T c = — / f(E)lm [h(E+) (A Vn (E + )) v ] dE, (25) 



where f(E) is the Fermi function, and E + = E + iO. 
Eq. (|25|), giving the Curie temperature as a function of 
the parameters of the system, is the main result of the 
present work. 

It is worth analyzing the limiting cases of this equation. 
In J ^> W limit, shifting the energy by J, we obtain Eq. 
@ in the form |E| 



T, 



Aw 

7T 



f(E)lm 



(9v n (E+)f v 



w(g Vn (E + ) 2 ), 



dE, (26) 



where gy n (E) — {E ~w (gv n )y — Vn)^ 1 . In the Appendix 
we compare the J <C W limit of Eq. ( p5f ) with the results 
of the Ruderman-Kittel-Kasuya-Yosida (RKKY) theory 



IV. PHASE DIAGRAM 



Consider first the the phase diagram (PD) of the sys- 
tem in case of no quenched disorder. The PD is presented 
on Fig. 1. 
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FIG. 1. The phase diagram for the case of no quenched 
disorder in the coordinates of relative strength of the Hund 
exchange J/W and electron concentration n. 



In the region, where Eq. (25), gives negative value 
for the T c , ferromagnetism is precluded at any tempera- 
ture. From our consideration we can say nothing about 
the nature of the non- ferromagnetic phase (or phases), 
but we know from the theory of the RKKY interaction 
]l6[ , that for small J/W (and no disorder), the ground 
state for the intermediate electron concentration is anti- 
ferromagnetic. One can say that the situation with finite 
Hund exchange is equivalent in some sense to the situa- 
tion with the infinite Hund exchange and antiferromag- 
netic superexchange |T^ ]. 

We consider the model of the disorder in which V n = V 
with the probability x, and V n = with the probabil- 
ity 1 — x, thus x being the concentration of impurities. 
Solving equation for the locator for the case of strong 
quenched disorder (V/W = 1 and x — .3) we obtain the 
PD, which is presented on Fig. 2. 
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FIG 2. The phase diagram for V/W = 1 and x = .3 in the 
coordinates of relative strength of the Hund exchange J/W 
and electron concentration n. 

It is interesting that ferromagnetism is now precluded 
in much larger region of the J/W — n plane. 

This research was supported by the Israeli Science 
Foundation administered by the Israel Academy of Sci- 
ences and Humanities. 



APPENDIX A: DMFA VS RKKY 

In J <C W limit from Eq. ( ^5| ) after simple algebra we 
obtain 



1 



—Im(g v (E + ) 2 ) v \dE 



(Al) 



where N V {E) = -(l/7r)Im gv(E+) is the DOS (per one 
spin direction) , and gy (E) is the solution of the equation 

1 



9v(E) = 



E-2w (g v ), 



V 



(A2) 
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In fact, in the weak exchange case we can expand Eq. 
(||) with respect to exchange interaction. The second or- 
der term, which is only important for us, for the diagonal 
matrix element of the Green's function G is 



where 



Gg>(IS)=E 



E-ei 



(A12) 



(2) 



G 



where 



N 



E 

I'nn' 



< ll'\nn' > 
(E - q) 2 (£ - ijO 



(m n • o-)(m„/ • a), (A3) 



(l/7rA0Im£„G 



(0) 



< ll'\nn' >= ^(n)iJi,(n)?Pp(n')iljL(n') 



(A4) 



and "0 are the eigenfunctions of the non-magnetic part 
of the Hamiltonian: (Ho + V)ipi — eitpi. The thermody- 
namic potential of the electron subsystem is given by the 
formula 



and NW(E) 

Eqs. (Al) and (All) look very much alike. The only 
difference between them is ensemble averaging in Eq. 
(Al) vs site averaging in Eq. 
averaging, that is NM(E) = 



flAUp . The DOS is self- 
{Nv{E)) v , because it in- 



rt 2 



1 



111 



-f}(E-^) 



(A5) 



Calculating Tr with respect to spin indices and integrat- 
ing by parts we obtain 



^ (2) = ^ £ In 



/Ulri ■ m 



where 



f(E)lm 



E 

. w 



volves the locator itself |t8|. Moreover, it is known that 
the CPA results for the DOS are reasonable even in the 
case of strong disorder. So for low electron concentration, 
when only the first term in Eq. (Al)(orEq. (All)) is im- 
portant, the equations arc equivalent. For higher electron 
Im Tr g[ ; 2 ' (E + )dE. concentration the term with the square of the locator de- 
i creases the T c . According to the RKKY theory, the Curie 

temperature goes through zero at approximately n = .25 
for the three principal cubic lattices |ll| (for the case of 
no quenched disorder). Eq. (Al) for this case gives crit- 
ical concentration n — .4. This comparison allows us to 
estimate the degree of agreement between the results of 
the DMFA and RKKY theory for the case considered. In 
the opposite case of very strong quenched disorder, the 
diffe rence between the square of the locator in Eqs. ( Al) 
and ( |All ) becomes even larger, due to the effects of local- 
ization, which are absent in the CPA. But the influence 
(A7)bf the localization on the destruction of ferromagnetism 
in the DE model demands additional consideration. 



< ll'\nn' > 



[E + -ei){E + -ev 



(A6) 



dE. 



The exchange integral, after fulfilling in Eq. (A7) inte- 
gration with respect to dE, can be presented as 

T t2 f( e l') ~ f( 6 l) - ;,/| / A o \ 

J-nri = J / < \nn > . (A8) 

w 1 1 

To see the connection between the RKKY approxima- 
tion and Eq. (|Al|), let us make in Eq. (A6) a mean field 
approximation m„ • m n > = M 2 . Thus obtained potential 
can be used to construct the Landau functional of the 



system, which gives |14 
Finally using the formula 



' 3 NM 2 ' 



(A9) 



£ < ll'\nn' >= 5 W - ]T \Mn)\ 2 \ipi'(n)\ 2 , (A10) 



we obtain 



2 J 



f(E) 



dN(°\E) 



dE 



dE, 



(AH) 
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